DESY 02-127 
September 2002 

Femtosecond Resolution Experiments at 
Third-Generation Light Sources: a Concept Based on the 
Statistical Properties of Synchrotron Radiation 

E.L. SaldinS, E.A. Schneidmiller S, M.V. YurkovS 

^ Deutsches Elektronen- Synchrotron (DESY), Notkestrasse 85, D-22607 Hamburg, 

Germany 

^ Joint Institute for Nuclear Research, Dubna, 141980 Moscow Region, Russia 
Abstract 

The paper describes a new concept of visible pump/X-ray probe/slow detector 
experiments that could be performed at third-generation synchrotron light sources. 
We propose a technique that would allow time resolution up to femtosecond capa- 
bilities to be recovered from a long (100 ps) X-ray probe pulse. The visible pump 
pulse must be as short as the desired time resolution. The principle of operation of 
the proposed pump-probe scheme is essentially based on the statistical properties 
of the synchrotron radiation. These properties are well known in statistical optics 
as properties of completely chaotic polarized light. Our technique utilizes the fact 
that, for any synchrotron light beam there exist some characteristic time (coher- 
ence time), which determines the time-scale of the random fluctuations. The typical 
coherence time of soft X-ray synchrotron light at the exit of monochromator is in 
the femtosecond range. An excited state is prepared with a pump pulse and then 
projected with a probe pulse onto a final ion state. The first statistical quantity 
of interest is the variance of the number of photoelectrons detected during syn- 
chrotron radiation pulse. The statistics of concern are defined over an ensemble of 
synchrotron radiation pulses. From a set of variances measured as a function of 
coherence time (inversely proportional to monochromator bandwidth) it is possible 
to reconstruct the femtosecond dynamical process. 
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1 Introduction 



Time-resolved experiments are used to monitor time-dependent phenomena. The 
study of dynamics in physical systems often requires time resolution beyond the 
picosecond capabilities presently available with synchrotron radiation. Femtosecond 
(10~^^s) capabilities have been available for many years at visible wavelengths. 
There exists a wide interest in the extension of femtosecond techniques to the soft 
X-ray, and X-ray regions of spectrum. 

The time resolution of pump-probe experiment is, obviously, determined by the 
duration of the pump as well as the resolution of the probing. The pump pulse must 
always be as short as the desired time resolution. In typical scheme of a pump- 
probe experiment (see Fig. ^ the short probe pulse follows the pump pulse at some 
specified delay. The signal recorded by the slow detector then reflects the state of 
the sample during the brief probing. The experiment must be repeated many times 
with different delays in order to reconstruct the femtosecond dynamical process 




The obvious temporal limitation of the typical visible pump/X-ray probe technique 
is the duration of the X-ray probe. Here we will concentrate on the performance 
obtained with third-generation storage rings. At these sources, the X-ray pulse du- 
ration is about 100 ps. The new principle of pump-probe techniques described below 
offers a way around this difficulty. We propose visible pump / X-ray probe technique 
that would allow time resolution up to femtosecond capabilities to be recovered from 
a long X-ray probe pulses (see Fig. ^). The principle of operation of the proposed 
pump-probe scheme is essentially based on the statistical properties of the syn- 
chrotron radiation. Synchrotron radiation is a stochastic object and at a given time 
it is impossible to predict the amount of energy which flows to a sample. For any 
synchrotron light beam there exist some characteristic time, which determines the 
time-scale of the random fluctuations. This characteristic time is called, in general, 
the coherence time Tc of the synchrotron light beam. Its magnitude is of the order 
of the inverse of the frequency spread of the beam. In all the theory that follows, 
attention is restricted to synchrotron light beams whose frequency spreads are small 
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Fig. 1. Typical scheme of a "short pump-short probe" experiment 
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Fig. 2. The new scheme of "pump-probe" experiment. The principle of operation of 
the proposed scheme is based essentially on the statistical properties of synchrotron 
radiation, which possesses all the features corresponding completely chaotic polar- 
ized light 

compared with the mean frequency, that is, where ujqTc is very much larger than 
unity. Figure ^ illustrates the type of fluctuations that occur in the synchrotron 
radiation beam intensity and phase. The figures have been constructed by a com- 
puter simulation of a synchrotron light source in which the summation is carried 
out explicitly for the real number of electrons {N ~ 10^'^) in the electron bunch. It 
is seen that substantial changes in intensity (and phase) occur over a time span Tc, 
but these quantities are reasonably constant over time intervals At <C Tq- 

Our technique utilizes the fact that, over the long probe pulse, time correlations 
exist in the probe pulse intensity. It was shown that the femtosecond timescale 
associated with intensity fluctuations in synchrotron light pulses make them well 
suited for time-resolved studies. The first statistical quantity of interest is the vari- 
ance of the detector counts. The statistics of concern are defined over an ensemble 
of synchrotron radiation pulses. From a set of variances measured as a function of 
correlation time it is possible to reconstruct the femtosecond dynamical process. 



2 The principle of pump-probe techniques based on the statistical 
properties of synchrotron light 



The shot noise in the electron beam causes fluctuations of the beam density which 
are random in time and space. As a result, the radiation, produced by such a beam, 
has random amplitudes and phases in time and space. These kinds of radiation 
fields can be described in terms of statistical optics, branch of optics that has been 
developed intensively during the last few decades and there exists a lot of experience 
and a theoretical basis for the description of fluctuating electromagnetic fields ||3|,^ . 

We will begin in this section by dealing with some general notions of statistical 
optics. Some of the statements will be quite precise, other only partially precise. 
In the framework of statistical optics the radiation field is characterized by notions 
such as time and space coherence. Let us illustrate these notions. To be specific, we 
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Fig. 3. Time-dependence of the intensity (upper plot) and phase (lower plot) for 
synchrotron light, obtained from a computer simulation. The dashed line shows the 
mean value of the intensity averaged over an ensemble of radiation pulses 

consider the radiation pulse of nearly monochromatic radiation having a duration 
and bandwidth equal to T and Auj, respectively. If T ^ 1/Alu, the radiation is 
partially coherent in time. The time coherence is of about ~ 1/Aui. The physical 
sense of this notion is as follows. Let us separate the radiation pulse by a splitter. 
These two pulses pass different path lengths and are then combined together. If 
the difference between the path lengths is less than cTc, we see the interference 
pattern at the end for a shot-to-shot averaging (averaging over a large number of 
radiation pulses). Using simple physical language, we can say that the radiation 
field is correlated within the time of coherence. 

The notion of space coherence can be explained in the same way. Let us direct the 
radiation onto a screen with two pinholes and look at the interference pattern in 
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the far diffraction zone when changing the distance between the pinholes. When 
the pinholes are located close to each other, we see a clear diffraction pattern. As 
the distance between the holes increases to some value, D > Arc, we come to 
the situation when we do not see an interference pattern after averaging over the 
ensemble of pulses. Qualitatively, the value of 7r(Arc)^ is referred as the area of 
coherence of the radiation pulse. The coherence volume is defined as the product of 
the coherence area, 7r(Arc)'^, with the coherence length, cTc- If one uses the notions 
of quantum mechanics, this volume corresponds to one cell in the phase space of the 
photons. The number of photons in the coherence volume is also referred to as the 
degeneracy parameter, 5^. Physically this parameter means the average number 
of photons which can interfere, or the number of photons in one quantum state 
(one mode). It should be noted that third-generation synchrotron light sources can 
produce soft X-ray radiation with very high degeneracy parameter and in such 
cases we can state that the classical approach is adequate for a description of the 
statistical properties of the synchrotron radiation. 

Electron beam is composed of large number of electrons, thus fluctuations always 
exist in the electron beam density due to the effect of shot noise. When the electron 
beam enters the undulator, the presence of the beam modulation at frequencies close 
to the resonance frequency initiates the process of radiation. The shot noise in the 
electron beam is a Gaussian random process. The monochromator can be considered 
as a linear filter which does not change the statistics of the signal. As a result, we can 
define general statistical properties of the radiation after monochromator without 
any calculations. For instance, the real and imaginary parts of the slowly varying 
complex amplitudes of the electric field of the electromagnetic wave have a Gaussian 
distribution, the instantaneous radiation power fluctuates in accordance with the 
negative exponential distribution, and the energy in the radiation pulse fluctuates 
in accordance with the gamma distribution. We can also state that the spectral 
density of the radiation energy and the first-order time correlation function should 
form a Fourier transform pair (this is so-called Wiener Khintchine theorem). Also, 
the higher-order correlation functions (time and spectral) should be expressed in 
terms of the first-order correlation functions. These properties are well known in 
statistical optics as properties of completely chaotic polarized radiation [^. This 
demonstrate one of the most beautiful things about statistical optics - how much 
can be deduced from so little. This brings up an interesting question: Why is it 
that the noise in the electron bunch is a Gaussian random process? We will give an 
explanation in section 3. 

The arrangement of components for a long probe pulse experiment is shown in Fig. 
In order to use the radiation from synchrotron light source it first has to be spatially 
filtered and monochromatized. We assume that synchrotron radiation reaching the 
sample surface has full transverse coherence. With this assumption, attention can 
be concentrated completely on temporal coherence effects. The coherence time of 
synchrotron light is inversely proportional to the monochromator bandwidth. So 
the narrower we make the monochromator exit slit, the wider the coherence time 
gets. Any measurement of photodetector counts will be accompanied by certain 
unavoidable fluctuations and we wish to determine the statistical distribution of 
the number of photocounts observed in any radiation pulse. 
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Fig. 4. Proposed pump-probe experimental setup 
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Fig. 5. Schematic ionization scheme for pump-probe ionization in the visible 
pump/X-ray probe transition. A typical experiment would be to couple two elec- 
tronic state |1) and |2) resonantly by a femtosecond pump- field. The dynamics is 
then probed by ionization and detection of the fluctuations of photoelectron counts 
per pulse as a function of the coherence time of the probe synchrotron radiation 
pulse 



Now we wish to consider one example which shows the relationship between the 
fluctuations of photocounts and the coherence time in a circumstances that easy to 
understand. We will not discuss applications of the proposed pump-probe scheme. 
There are, in fact, many applications that could be discussed. We have chosen for 
emphasis here experiments aimed at measuring the lifetime of electronic excitation 
in condensed matter. For our present purposes we would like to imagine a somewhat 
idealized experiment. Several assumptions are made about the process of internal 
conversion of the electronic excitation. First, the time dependence of the electronic 
occupation is assumed to be known a priory - this is the exponential dependence 
nit) = TT-o exp(— t/ro). In practice, this is a good assumption. The purpose of the 
measurement is presumably to determine the lifetime tq. The typical lifetime of 
the electronic excitation is in the femtosecond range. This type of time-resolved 
experiments may be understood as preparing an excited state with a pump pulse and 
then projecting it with a probe pulse onto a final state. The final state should be well 
characterized. Photoionization detection has several advantages. The ion state may 
be well characterized by independent methods such as high resolution photoelectron 
spectroscopy. Schematic ionization scheme for pump-probe ionization is illustrated 
in Fig. ^. A typical experiment would be to couple two electronic state | 1) and 
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I 2) resonantly by a femtosecond pump-field at a frequency uj in the visible range. 
The dynamics is then probed by ionization and detection of the photoelectrons 
as a function of the coherence time of the probe synchrotron radiation pulse. A 
time-resolved experiment inherently begins with the initiation of the process under 
study, at some more or less accurately defined instant in time. We assume that the 
pump pulse can be approximated by a (5-function. With that simplifying assumption, 
the number of photoelectrons K detected during a synchrotron radiation pulse is 
directly proportional to the integrated value of the physical parameter (occupation 
number) to be measured: 



K (X j I{t)no exp(-t/ro) d t 



where I{t) is the instantaneous intensity of synchrotron radiation. At time t = 
the pump pulse perturbs the sample. Here we assume that the pump-probe process 
under consideration can be separated from what might be called parasitic processes. 
In practice, not every count will be detected, and some false counts will register. We 
are not considering such details because our interest is in the fundamental aspects 
of the problem. 

A quantity of considerable physical interest is the variance of the photocount dis- 
tribution 



Determination of the variance cr^ as a function of the monochromator bandwidth 
gives us information on the lifetime tq. Why? First, for a lifetime tq that is 
much shorter than the coherence time Tc of the probe pulse, the number of photo- 
electrons is to an excellent approximation, simply the product of the instantaneous 
intensity and the lifetime tq, 

oo 

K oc I{t)no exp(— t/ro) dt ^ /(O)noro for AwmTo < 1 • 



Within the scaling factor, therefore, the probability density function of K is approx- 
imately the same as the density function of the instantaneous intensity. Using the 
well-known results obtained in the framework of statistical optics, we can state that 
when the lifetime is much shorter than the inverse of the monochromator bandwidth, 
the variance tends to unity: 

(Tk ~ 1 for AwmTo < 1 . 

At the opposite extreme, with a lifetime much longer than the coherence time, the 
fact that many independent fluctuations of the instantaneous intensity occur within 
the interval tq implies, according to the theory of completely chaotic light, that the 
variance is inversely proportional to the monochromator bandwidth. So, we obtain 
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Fig. 6. Typical behavior of the variance associated with the measurement of the 
number of photoelectrons per pulse as a function of lifetime and monochromator 
bandwidth 

the rule that the product of lifetime and monochromator bandwidth is of the order 
of 

Auj^tq ~ fjj^^ for AWniTO > 1 ■ 

This relation is illustrated in Fig. ^, which shows the parameter fij^^ plotted against 
AwmTQ. The exact results can in fact be found for certain monochromator line 
shapes using the formalism introduced in Section 2. The dependence of on the 
exact shape of the signal and of the monochromator profile is rather weak and can 
be ignored outside the range 0.1 < Auj^tq < 10. 

It is important to note that the photocount fluctuations discussed above originate 
from the shot noise in the electron beam. Nevertheless, the problem of photocount 
fluctuations is more complicated and there is another (parasitic) fluctuation effect 
that can be important for the proposed pump-probe experiments. The problem is 
the time jitter of electron bunch and pump pulse. Any measurement of photocounts 
will be accompanied by these unavoidable fluctuations. To solve this problem, fem- 
tosecond pump pulses from the laser system should be synchronized to the master 
clock of the storage ring rf system with phase-locking technique. At present the 
best achievable synchronization between the optical laser and the synchrotron light 
source is of the order of a few picoseconds. On the other hand the electron bunch 
length is about of 100 ps. Assuming that the electron bunch envelope is Gaussian 
distribution (in practice, this is a good assumption) we can find that the parasitic 
photocount variance reaches a value of about a per cent only. Thus we find that 
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pump pulse jitter should not be a serious limitation in our case. 



The next problem is to define the sensitivity of the proposed pump-probe method. 
A general study of the noise fluctuations associated with the photoevents would be 
nontrivial. The difficulty arises in simultaneously including the effects of both the 
classically induced fluctuations and Poisson noise fluctuations associated with the 
basic interaction of light and matter. The ratio of the classical variance to the "pho- 
ton noise" variance is named as a photocount degeneracy parameter 6c (see section 
4). Physically speaking, the count degeneracy parameter can be interpreted as the 
average number of counts that occur in a single coherence interval of the incident 
(transversely coherent) radiation. It can also be described as the average number of 
counts per "degree of freedom" or per "mode" of the incident wave. When Sc <^ 1, 
it is highly probable that there will be no more than one count per coherence inter- 
val of the wave, with the result that quantum noise predominates over classically 
induced noise. On the other hand, when 5c ^ 1, there are many photoevents pre- 
sented in each coherence interval of the wave. The result is a "bunching" of the 
photoevents by the classical intensity fluctuations, and an increase of the variance 
of the counts to the point where the classically induced fluctuations are far stronger 
than the quantum noise variations. 

Let us consider the case when the degeneracy parameter is much larger than unity, 
(5c ^ 1- The fluctuations of the photoevents are defined mainly by the classical noise 
in this case and the signal-to-noise ratio for the proposed device can be written in 
the form 

N s/{{K - (K)}*) - {(K - (K))-')-' ^' 

where A'shot is the number of independent measurements averaged in the accumu- 
lator (total number of shots). The determination of the signal-to- noise ratio is a 
nontrivial problem. Fortunately, in the particular case, namely, Tc ^ tq, a much 
simplified analysis is sufficient. In this case the expression (|l|) can be reduced to 



S/N ~ ^A^shot/S for (5c > 1 , Aco^to < 1 . 



The signal-to-noise ratio depends on the degeneracy parameter when 5c is much less 
than unity. In this case one can derive that 

S/N ~ S^J^^N^t for (5c < 1 , Aco^To <C 1 . 

The significance of the result we have obtained cannot be fully appreciated until 
we determine typical values of the degeneracy parameter that can be expected in 
practice. Let us consider the case of radiation with a sufficiently long coherence time 
that Tc ^ tq. Under such condition, the photocount degeneracy parameter can be 
estimated simply as (see section 5) 

5c ^ r?(AWmTo) |-Rm '^ ^4^ 1 ^^"^ ^^mTo < 1 , (2) 
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where r] is the quantum efficiency of photoelectron production, is the monochro- 
mator throughput, A is the radiation wavelength, -Bpeak is the peak spectral bright- 
ness. Our physical interpretation of the result (|2|) is as follows. In order to use the ra- 
diation from the undulator it first has to be spatially filtered and monochromatized 
(see Fig. Q). We interpret the factor R^X^ Bpes,k/ {4c) as representing the number 
of photons per coherence interval at the monochromator exit, whereas the factor 
AcUmTo represents a ratio of integration time tq to coherence interval Tc ~ l/Awm- 
Let us present a specific numerical example for the case of a third-generation syn- 
chrotron light source. The spectral brightness delivered by an undulator s at a wave- 
length of A ~ 10 nm is about Speak - 10^^ photons/s/0.1%BW/ mm^/mrad^ and 
the degeneracy parameter is about 6c — 10^?7i?m for the case when Tc ~ tq. 

As we mentioned above, the dependence of o"^^(tc) on the exact shape of the signal is 
rather weak. As a result, the behaviour of the variance as a function of synchrotron 
radiation coherence time provides only very little information on the sample dy- 
namics. Nevertheless, the described pump-probe technique can be used to quote a 
"characteristic time" of the sample dynamics. The typical procedure is to assume a 
signal shape (generally exponential shape), and to "determine" the relaxation time 
Tq from the known ratio between the tq/tc and cj^. 

A more complex instrument can be used to characterize the shape of the signal 
function. The pump-probe technique based on a correlation principle we are describe 
in section 7 is to use two spatially separated samples and record correlation of the 
count fluctuations for each value of the delay r between two pump laser pulses. 
Averaged (shot to shot) count product contains information about the shape of the 
signal function. Another advantage of the correlator measurement is the possibility 
to remove the monochromator between spatial filter and sample and thus to increase 
the count degeneracy parameter. 



3 Statistical properties of synchrotron radiation 



The present section considers the statistical properties of the intensity fluctuations 
in synchrotron light. In order to give the subject a semblance of continuity, it will be 
desirable to introduce considerable matter which will be found in any of standard 
works on statistical optics. Much of the credit for stimulating modern development 
in statistical analysis of synchrotron radiation is due to Professor's Goodman book 
[Q, which laid the groundwork for the principles of pump-probe techniques based 
on statistical properties of synchrotron light. Quite basic similarities between the 
methods described here and the methods introduced into X-ray free electron laser 
physics by authors in ||5| should be recognized. We here consider the classical the- 
ory of the fluctuation experiment, reserving the quantum discussion for section 4, 
where the theory of intensity fluctuations is reconsidered in terms of the quantized 
radiation field. 
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3.1 Shot noise in the electron beam 



In this subsection we study the statistical properties of the shot noise in the electron 
beam. It should be noted that the process under study is nonstationary with finite 
pulse duration, so in what follows the averaging symbol (• • •) means the ensemble 
average over bunches. 

Let us consider the microscopic picture of the electron beam current at the en- 
trance of the undulator. The electron beam current is made up of moving electrons 
randomly arriving at the entrance to the undulator: 

N 

J{t) = {-e)Y^5{t-tk) , 

k=l 

where 5{- ■ ■) is the delta function, (— e) is the charge of the electron, is the number 
of electrons in a bunch, and is the random arrival time of the electron at the 
undulator entrance. The electron bunch profile is described by the profile function 
F(t). The beam current averaged over an ensemble of bunches can be written in the 
form: 

(J(t)) = {-e)NF{t) . 

For instance for an electron beam with Gaussian distribution of the current along 
the beam, the profile function is: 



F{t) 



The probability of arrival of an electron during the time interval {t,t + d t) is equal 
to F{t)dt. 

The electron beam current, J(t), and its Fourier transform, J{uj), are connected by 

oo jy 

= I e'^'Jit) d t = (-e) J2 , (3) 




°o ^ 



k=l 

k=l 



It follows from these expressions that the Fourier transform of the input current, 
J{lv), is the sum of a large number of complex phasors with random phases 4>k = ^tk- 
When the electron pulse duration, ctt, is long, ojai » 1, the phases (f)^ can be 
regarded as uniformly distributed on the interval (0, 2ti). The formal summation of 
the phasors is illustrated in Fig. |^ The calculation is an example of the "random 
walk" problem, well known in the theory of stochastic process. In this case we can 
use the central limit theorem and conclude that the real part and the imaginary part 
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Fig. 7. Diagram to show the amphtude and phase of the resultant vector (total 
complex amplitude) formed by a large number of unit vectors (elementary com- 
plex-valued contributions exp(ia;tfc)), each of which has a randomly chosen phase 
angle 

of J {to) are distributed in accordance with the Gaussian law The probability 
density distribution of | J{oj) p is given by the negative exponential distribution: 



1 



(I J{^) 



exp 



(I J{^) 



(4) 



3.2 First-order spectral correlation 



Let us calculate the first-order correlation of the complex Fourier harmonics J(w) 
and J{uj'): 

N N 

{J{u)r{J)) = e^{J2J2 exp(iu;tfc -iio'tn)) ■ 

fc=ln=l 

Expanding this relation, we can write: 

N 

{J{lj)J*{J)) = e^(^ exp[i(u; - cj')*fc]) + X!^*^^P(i'^^fc))^^^P(~i^'*")) ' 

fc=l k^n 

Expression (exp(ia;tfc)) is equal to the Fourier transformation of the bunch profile 
function F{t): 



(exp(ia;tfc)) = / F(tk)e''^'^ dtfc = F{uj) . 
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Thus we can write: 



{J{uj)J*{uj')) = e^NF{uj - J) + e'^N{N - 1)F{uj)F{uj') . 
When 

N I F(cj) I V 1 , 

we can write the following expression for the first-order spectral correlation: 

{J{uj)J*{uj')) =e^NF{u-uj') . 
The Fourier transform of the Gaussian profile function has the form: 



(5) 



(6) 



F{u)) = exp 



2^2 



CO a. 



For the specific cases of a Gaussian profile of the electron bunch, the first-order 
correlation of the complex Fourier harmonics, J{uj) and J {to'), has the form: 



{J{uj)J*{uj')) = e^iVexp 



(uj - Lo')'^a^ 



(7) 



Let us discuss the region of validity of the approximation (|5|) . The physical meaning 
of @ is that the frequency oj has to be large enough, ooa^ S> 1. To be specific, we 
consider a numerical example for a Gaussian bunch profile. The value of | F{uj) ^ is 
equal to exp(— 10^^) when a^uj = 10®. So condition (^) is always fulfilled in practice. 



3.3 Second-order spectral correlation 



Let us calculate the second order correlation of the complex Fourier harmonics J{io) 
and J{lo'): 

N N N N 

(I JH |2| JV) |2) = e^^ ^ ^^exp[i.;(i„-U+ia;%-g]) . 

n=l m=l p=l q=l 

The A^^ terms in this sum can be set in 15 classes Q. When condition (^) is fulfilled, 
only two of them are of importance, corresponding to (n = m, p = q, n ^ p) 
and {n = q, m = p, n ^ m). Thus we can write: 

(I J» PI J{J) |2) = (I J» J^J) |2)+ I (J»J(a;')> ? ■ (8) 

Substituting (0) into (|8|) we obtain: 

(I J{uj) Pi J{J) |2) = e^Ar2(i+ | f{u - J) p) . (9) 
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3-4 The origin of shot noise in the electron beam 

Why is it that the noise in electron beam is a Gaussian random process? Why is 
this the right rule, what is the fundamental reason for it, and how is it connected to 
anything else? The explanation is deep down in quantum mechanics S. It is gener- 
ally accepted that quantum mechanics provides the best current picture of physical 
phenomena, and the most complete description of the radiation field must be sought 
in quantum-mechanical terms. The quantum theory of radiation predicts the exis- 
tence of zero-point electromagnetic field. The first step is to show that a field mode 
is equivalent to harmonic oscillator. To describe a field mode quantum mechani- 
cally, we simply describe the equivalent harmonic oscillator quantum mechanically. 
The vacuum state | 0) has no photons, but predicts that it nevertheless has an en- 
ergy hLo/2. From the probability distribution for a ground-state harmonic oscillator, 
we easily obtain the probability distribution for the electric field in vacuum state 
I 0). What is the physical significance of these vacuum-state expectation values? 
One thing they indicate is that the electromagnetic vacuum is a stationary state of 
the field with statistical fluctuations of the electric and magnetic fields. As far as 
measurements are concerned, however, it is often argued that the entire universe is 
evidently bathed in zero-point electromagnetic field, which is distributed according 
to a Gaussian probability distribution 

Relativistic electron beam passing through bending magnets emits synchrotron ra- 
diation, a process that leads to damping. In standard classical electrodynamics there 
is only the radiation reaction field to act on a single particle in the vacuum. All six 
degrees of freedom for electron motion in storage ring are damped and the actual 
value of electron beam emittance would be equal to zero. We deduce the synchrotron 
radiation formula using a classical treatment. So, the reasonable question arises of 
whether such a classical model describes correctly the synchrotron radiation in the 
storage ring. Simple physical considerations show that such a description is valid for 
any practical situation. The classical electrodynamics approach can be used when 
the energy of the radiated photon huj is much less than the energy of the electron. As 
a rule , fioj / {mc?^) ~ 10^^ and quantum effects are negligible. Classically we gen- 
erally assume implicitly that the homogeneous solution of the Maxwell equations is 
that in which the electric and magnetic fields vanish identically. That is, we assume 
that there are no fields in the absence of any sources. In the absence of sources the 
vacuum field is simply zero, with no energy or fluctuations whatsoever. This is not 
to say that we cannot have source-free fields in classical electrodynamics. Rather, 
the absence of zero-point fields in standard classical theory lies in the assumption 
that there are no fields in the absence of any sources. 

We can go beyond standard classical theory and postulate the existence of zero-point 
electric and magnetic fields in the absence of any sources. The difference between 



This subsection is rather abstract side tour. Reader can therefore, skip over it 
and continue with subsection 3.5. Noise in electron beam is a Gaussian random 
process. For the moment, reader will just have to take it as one of the rules of the 
world . 
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classical and quantum electrodynamics, together with the evident importance of the 
fluctuating vacuum field in quantum electrodynamics, suggests the adoption of a 
different boundary conditions in classical electrodynamics: instead of assuming that 
the classical field vanishes in the absence of sources, we can assume that there is 
a fluctuating classical field with zero energy hoj/2 per mode. As long as this field 
satisfies the Maxwell's equations there is no a priori inconsistency in this assump- 
tion. The appearance of fi in this modification of classical electrodynamics implies 
no deviation from conventional classical ideas, for fi is regarded as nothing more 
than a number chosen to obtain consistency of the predictions of the theory with 
experiment. The resulting theory is able to account quantum effects within a fully 
classical framework [||. 

We note that electron storage rings have provided strong evidence for the reality of 
zero-point energy. In the experiments with stored electrons the zero-point energy 
sets a lower limit to freezing of the electron motion. If there were no noise associated 
with the vacuum fields, the actual value of electron beam emittance would be equal 
to zero. In storage rings damping will counteract quantum excitation leading to an 
equilibrium. If one looks over the derivation of the quantum diffusion of the electrons 
in the storage ring, it becomes clear that the argument can be couched in the 
language of stochastic electrodynamics rather than quantum electrodynamics. That 
is , all that is really required in that derivation is zero-point energy fiLo/2 per mode of 
the electromagnetic field. Whether this zero-point energy is of quantum or classical 
origin is irrelevant for the purpose of deriving the diffusion coefficient - stochastic 
electrodynamics accounts perfectly well for the quantum diffusion of electrons in 
storage ring. In this case both the field and the particles are treated classically and 
we have no practical need in accelerator physics for quantum electrodynamics. 

The diffusion terms will introduce a statistical mixing of the phases and after some 
damping times any initial azimuthal variation of the phase space density will be 
washed out. Electrons in a storage ring are evidently bathed in zero-point electro- 
magnetic field, which is distributed according to a Gaussian probability distribution. 
This is the explanation of the relation between vacuum fluctuations and shot noise 
in the electron beam. 



3.5 Analysis of synchrotron radiation properties in the frequency domain 

Above we described the properties of the input shot noise signal in the frequency 
domain. The next step is the derivation of the spectral function connecting the 
Fourier amplitudes of the output field and the Fourier amplitudes of the input 
noise signal. In the first analysis of the problem, we adopt some rather simplifying 
assumptions that are only occasionally met in practice. Following this simplified 
analysis, however, we show how the validity of the results can be extended to a far 
wider range of conditions than might have been thought at the start. 

We will investigate the synchrotron radiation in the framework of the one-dimensional 
model. For simplicity we consider linear polarization of synchrotron radiation, and 
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use a scalar representation of the radiation field. Nevertheless, all the results are 
valid for any polarization. The one-dimensional model describes the radiation of the 
plane electromagnetic wave 

Ey{t) = Ey{t) eyip[iuj{z/c - t)] + C.C. , 

by the electron beam in the planar undulator. The electric field of the electromag- 
netic wave in the time domain, Ey{t), and its Fourier transform, E{uj), are connected 
by 

oo 



When u < the Fourier harmonic is defined by the relation E*{lo) = E[—lo). The 
Fourier harmonic of the transversely coherent electromagnetic field at the spatial 
filter exit and the Fourier harmonic of the current at the undulator entrance are 
connected by the relation: 

E{iv) = A{u:)J{lo) , tj > , (10) 

where A{ijo) is the spectral function of the undulator. For an undulator of periods 
each electron oscillates through A'„ cycles of its motion and thus radiates a wavetrain 
consisting of A'„ cycles of the electric field. The Fourier transform of this waveform 
which gives the spectral content of the fields, is sin(x)/x, where x = N^t[{/\uj /ujq) 
and Aw = [u) — loq) is the frequency shift away from the central maximum at uiq. 
Thus we can write (see section 4): 

2 . sin^(jVw7rALj/a;o) 

It is relevant to make some remarks on the region of applicability of the one- 
dimensional theory. One-dimensional model assumes the input shot noise and out- 
put radiation to have full transverse coherence. This assumption allow us to assume 
that the input shot noise signal is defined by the value of beam current (|^). In 
reality the fluctuations of the electron beam current density are uncorrelated in the 
transverse dimension. Using the notion of the beam radiation modes, we can say 
that many transverse radiation modes are radiated when the electron beam enters 
the undulator. The one-dimensional model can be used for the calculations of sta- 
tistical properties of transversely coherent synchrotron radiation. In practice such 
an assumption is valid for synchrotron light at the exit of a spatial fllter. During the 
spatial filtering process, the number of transverse modes decreases, and the contri- 
bution of the coherent radiation into the total radiation power is increased up to 
full coherence. With this assumption, attention can be concentrated completely on 
temporal coherence effects. 

Here we have used the simplest model of synchrotron light source with zero energy 
spread into electron beam. In addition, it is assumed that the electrons move along 
constrained sinusoidal trajectories in parallel with the undulator axis. The first 
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assumption is primarily a statement that the energy spread, A7/7, is so small that 
N^Aj/'y <C 1. Such a condition is generally valid in practice. The second assumption 
is not generally valid in practice, but it will be removed in section 5. 

Above we studied the properties of the Fourier harmonics of the shot noise in the 
frequency domain. The Fourier harmonics of the output radiation field are connected 
with the Fourier harmonics of the input shot noise by (|10|). It follows from 1^ that 
statistical properties of the Fourier amplitudes E{lli) are defined by the statistical 
properties of the Fourier amplitudes of the input current J{lo). In particular, it 
follows immediately from (^) that | E{u>) p is distributed in accordance with the 
negative exponential probability density function: 

pi\ m ?) = iTT^ -p . (11) 



For many practical applications of the synchrotron radiation a monochromator has 
to be installed at the undulator exit. We denote the frequency profile function of the 
monochromator as Gjn{uj). The linearity of Maxwell's equations and the fact that 
a monochromator can be treated as a linear filter allows one to write the Fourier 
components of the electric field of the synchrotron radiation in the following form: 

E{u;) = A{u;)G^{uj)J{uj) , uj > , (12) 

Let us calculate the correlation of the complex Fourier harmonics E{uj) and E{uj'). 
When the monochromator bandwidth is much less than the spectral bandwidth of 
the undulator radiation we obtain: 

{E{u;)E*icj')) =1 A{uo) P G,^iu - u;o)G*^iu;' - LOa){J{u:)r{u:')) . (13) 
It is seen from (|6|) that one can write (^) as 

{E{u)E*{uj')) = e^N I A{uJo) p G^{uj - ujq)G*^{uj' - ojq)F{oj - J) . (14) 
The first-order spectral correlation function is defined as 

{E{u)E*{J)) 

gi{uj-uj) = (15) 
^{\E{u:)\^){\E{u:')\^) 

Substituting (||) into (|l5|), we obtain 

I gi{LO-J) 1 = 1 F{lo-J) I . (16) 

The explicit expression for the first-order spectral correlation function of the syn- 
chrotron radiation emitted by the electron bunch with Gaussian profile has the 
form: 



gi{LO - u') = exp 



(17) 
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We define the interval of spectral cofierence as 

oo 

Awc = J I gi{Au) 1 2 d(Ac^) . 



(18) 



The value of the spectral coherence for the synchrotron source with a Gaussian 
electron bunch is given by 



TT 



F{Auj) 1^ d{Auj) = ^ 



(19) 



The second-order spectral correlation function is defined as 
(I Ei^) PI Eico') |2) 



52 (W - io') 



(I E{u;)\'){\ E{co')\^) 



(20) 



Using (^), (|l2D and (20) we obtain that the first and second order correlation func- 
tions are connected by the relation: 



g2{uj - u') = 1+ \ gi{uj - uj') 



r\ |2 



(21) 



which is also a general property of completely chaotic polarized radiation. The 
explicit expression for the second-order spectral correlation function for the syn- 
chrotron source with Gaussian electron bunch has the form: 



52 (w - cu') = 1+ I F(uj -Lo')\=l + exp -{uj - Jya^ 



3. 6 Energy fluctuation experiment 



In this subsection we are going to discuss the application of statistical optics to 
a practical device. We will find that many of the features of this specific problem 
are quite common in the general statistical theory of synchrotron light, and we will 
learn a great deal by considering this one problem in detail. 

Figure ^ shows the fluctuation experiment geometry under consideration. The prob- 
lem is a description of the fluctuations of the energy of the radiation pulse W at 
the detector installed after the spatial fllter and monochromator. The experiment 
is here analyzed in some detail to determine the conditions under which the fleld 
fluctuations, resulting from the chaotic nature of the light source, affect the energy 
fluctuations. Note that in this experiment, both monochromator resolution and ra- 
diation pulse duration play a role. There is no X-ray monochromator providing the 
coherence time Tc — A^/cAA comparable with radiation pulse duration (100 ps). 
The model experiment ignores these technical limitations. 

Using the expression for Poynting's vector and Parseval's theorem, we calculate the 
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spatial monochromator detector 
filter 

Fig. 8. Set-up of model experiment 
radiation energy in one radiation pulse: 



W 



CO oo 



where S is the transverse area of the detector. The energy, averaged over an ensem- 
ble, is given by the expression: 

oo 

W = 4^/(1 ^HP)d^, (22) 



where (| E{ijj) p) is calculated using @ and (|l3|): 

(I E{u) |2) =1 A{lo) |2| G^{uo) |2 (I |2) = e^N \ A{u;) \^\ GM I' • (23) 



Using (22) and (p3D, we calculate the average energy measured after the monochro- 
mator: 

oo ^ oo 



It is seen that the average energy is a function of the frequency profile of the 
monochromator and undulator. 

The variance of the energy distribution is calculated as follows: 



oo oo 



W /T/l/\2 oo _ oo _ 

Jdu:{\E{c.)\^)Jdu;'{\E{u')\^) 



Using definition (|20| ) of the second-order correlation function and (^), we reduce 
this expression to 

oo oo _ _ 

/d^/d^'(|^(a;)P)(| E{u:')\^)\gi{u:-uj')\' 

c^w = 35 : 5s : • (24) 



The analysis of this expression shows that the variance of the energy distribution 
after the monochromator is a function of the frequency profile of the monochroma- 
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tor, of the frequency profile of the undulator, and of the electron bunch formfactor 
F{uj-uj'). 



Let us consider the case of an electron bunch with a Gaussian profile and a monochro- 
mator with a rectangular line: 

I Lr^[UJ — UJo) 1=1 for \ LO — UJQ \ < 



G„i{oj — uJo) p= for \ uj — ujQ \ > 



2 



For simplicity we assume that the bandwidth of the monochromator is small and 
A{uj) is constant within the monochromator bandwidth. Then the integration of the 
expression (p^) provides the following result: 



= IT— erf — r— - - -r— O - exp 



-IT 



\ AuJr 



where erf(x) is the error function and Awc is given by expression (19). 



(25) 



Let us study the asymptotic behavior of ( |25D . When the monochromator bandwidth 
is much less than the interval of spectral coherence, the normalized dispersion tends 
to unity: 

~ 1 for Atjjn ^ AtJc • 

When the monochromator bandwidth is much larger than the interval of spectral 
coherence, the variance is inversely proportional to the monochromator bandwidth: 

2 AWc 

The next practical problem is to find the probability density distribution of the 
radiation energy after the monochromator, p(W). Using the well-known results ob- 
tained in the framework of statistical optics, we can state that the distribution of the 
radiation energy after the monochromator is described rather well by the gamma 
probability density function: 

where T{M) is the gamma function of argument M, and 1/M = cj^. This distribu- 
tion provides correct values for the mean value of W and for the variance ■ 



w 



I wp{w) dw = {w) , I J^^^' piW) dW = ^ 



u u 

The parameter M can be interpreted as the average number of degrees of freedom 
(or modes) in the radiation pulse. It follows from (|2l) that this parameter cannot be 
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less than unity. When M tends to the unity, (2(:) tends to the negative exponential 
distribution (pi]). For large values of M the distribution ( p6D tends to a Gaussian 
distribution. 



3. 7 Analysis of the radiation properties in the time domain 



An expression for the electric field of the electromagnetic wave as a function of time 
t can be obtained using the Fourier transform of It is convenient to isolate 
explicitly the slowly varying complex amplitude: 



Ey{t) = Eit)e-'''°^ + C.C. , (27) 



where 



oo 







Let us consider the specific case of a Gaussian monochromator line profile 



G'm(t^) I = exp 



2al 



(28) 



allowing calculations of the complex amplitude E{t) of the electric field at the 
monochromator exit in the limit of small monochromator bandwidth. In other 
words, we assume A{uj) to be constant within the monochromator bandwidth. 



Using (|IOD, (I), (|2§) and (0), we obtain: 

^ ^ (29) 



f 1 ^ 1 

E{t) = Eo <^—Yl^exp{iuJotk) exp[-al,{t - tk)'^] | . 



Since Um <C ujq, we can approximately let the amplitude and the phase of each 
random phasor contributing to the sum in ( |29| ) be independent of each other, and 
the phases be uniformly distributed on the interval (0, 2tt). Thus, the distribution of 
the instantaneous radiation power density ( oc | E{t) p) is the negative exponential 
distribution: 

p(| E{t) = -J^ exp (- ' ] . (30) 

The correlation between the radiation fields at times t and t' has the form: 

{E{t)E*{t'))eM-^Mt-t')] 

oo oo 

= ^y"du; J duj' exp{- itot + iu;'t'){E{uj)E*{u;')) 





20 



oo oo 



^ 4~2 J J du;' exp{-iuit + iuj't') 



x[(|^M|2)(|^(^')P)]'/'i^(^-'^')- (31) 



Here ( |15D and ( pi3[) have been used when rewriting the integral. In the present consid- 
eration we use the approximation which assumes the interval of spectral coherence 
is much less than the monochromator bandwidth 

AlVc <C dm . (32) 

Thus we can simplify the integral in the following way. We replace the expression 
in square brackets by (| E{uj) p), and after integration over to and Aa; = {to — uj') 
we obtain: 



{E{t)E*{t'))eM--^Mt-t')] 

oo oo 

- —j d(Au;)exp(-iAa;t)F(ALj) ^ dw exp[- ia;(t - t')](l ^('^) 1^) 



47r2 

m 

27r 







y"dcuexp[-iu;(t-t')](l E{uj) \^) , (33) 





where F{t) is the radiation pulse profile. We define the first-order time correlation 
function as follows: 



Using (|33[), we can write 

oo _ 

/ dw(| E{uj) |2) exp[- i(w - UQ){t - t')] 

9l{t - t') = ^ 5S 

/du;(|^(^)P) 



Since we deal with a narrow-band signal, the latter expression may be rewritten as 
follows 

oo 

/ d(Au;)(| E[Auj) |2)exp[-i(Ac^)(t-0] 

gi{t-t') = — 55 , (34) 

/ d(Aa;)(| ^(A^) P) 

— oo 

where Aa; = (w — loq). Therefore, the slowly varying correlation function and the 
normalized spectrum of the narrow-band signal are a Fourier transform pair. Re- 
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membering relation (0), we rewrite (^) in the following way: 

oo 

/ d{Au) I G^iAu) |2 exp[-i(Aw)(t - t')] 

9i{t - t') = --^ . 

/ d(ALu) I G^(Aa;) P 

It is seen from this expression that the first-order correlation function possesses the 
property gi{t — t') = g\{t' — t). When the monochromator line profile is symmetrical 
the function gi is real. For instance, for a monochromator with Gaussian line profile 
the correlation function gi is 

5i(r)=exp(^-^^ , (35) 

where r = (t — t'). Following the approach of Mandel, we define the coherence time, 
Tc, as 



oo 



Tc = / I giir) f dr . 



-oo 



For the case of the Gaussian line profile , the explicit expression for the coherence 
time is 



Tc = ^ . (36) 

(Tin 

The correlation between the radiation intensities at times t and t' is defined as 

(I m PI E{t') p) 

oo oo oo oo 

= YEtt^ y"dwi J duJ2 J duJ3 J d^^4exp[-i{uJi-uJ2)t-i{u;3-uj4)t')] 


x{E{u;i)E{oJs)E*{oj2)E*{u;i)) . (37) 

Taking into account (0), (^) and (^, we simplify the correlation in the integrand 
of (^) in the following way: 

= {E{ui)E*{uj2)){E{uJ:i)E*{uj^)) 

+{E{u;i)E*{u;i)){E{i03)E*{u;2)) , (38) 



and present the integral ( |37D as the sum of two terms. To calculate the integral, one 
should take into account ([Tsl), (|l6|), (|33|) and (|3^) , which leads to 



(I Eit) PI Eit') p) = (I E{t) p)(| Eit') p)+ I {Eit)E*it')) p . (39) 
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The second-order time correlation function is defined as follows: 



jjemiioPL, (40) 

(I E(t) |2)(| E(t') P) 

It follows from (p9|) and definitions of gi{t — t') and g2{t — t') that 

g2{t-t') = l+\gi{t-t')\^ . (41) 

The analysis of the obtained relations (^) and (pl|) shows that the synchrotron 
radiation possesses all the features of completely chaotic polarized light (see also 



(11) and (|2l|) obtained in the frequency domain). 



3.8 Statistical properties of the integrated intensity 



In a variety of problems, related with photon-counting statistics, finite-time integrals 
of instantaneous intensity occur. In many cases the fluctuations in the intensity 
are too rapid for direct observation, and what is measured is some average of the 
fluctuations over detector response time. The calculation can be performed without 
great difficulty in two limiting cases, namely, the response time very long and very 
short compared with the electron pulse duration. In section 3.6 we discussed in 
detail fluctuations of the energy of the radiation pulse at the slow detector installed 
after the monochromator (see Fig. ^). Note that for this class of experiments we 
actually assume that the integration time is much longer than the electron pulse 
duration. 

Next let us consider the opposite extreme. For our simplified discussion here let us 
assume that the integration time is much shorter than the electron pulse duration, 
leaving the electron pulse envelope unaffected. There is no photodetector providing 
a time resolution much shorter than the X-ray pulse duration. Present-day X-ray de- 
tection will be limited in time resolution of about 100 ps. The situation is somewhat 
different in pump-probe experiments. A pump-probe experiment has analogies with 
a femtosecond photodetector measurement. The basic analogy is that a femtosec- 
ond pump pulse interacts with the sample and excites it into non-equilibrium state. 
The sample thereafter relaxes towards a new equilibrium state. Typical relaxation 
times are in the femtosecond range. This process can be mapped by sending a X- 
ray probe pulse onto the sample. Such a device could be discussed as a femtosecond 
X-ray detector by recording a secondary process such as photoelectrons. 

Several assumptions are made about the profile of the photodetector gating function. 
For the present, consideration is restricted to the simplest case from the analytical 
point of view, namely, completely chaotic polarized radiation and an ideal detector 
with a stepped profile gating function. So, the next problem is the description of 
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the fluctuations of the radiation energy W detected during a finite time interval 5T: 

t+5T 



w= J I{t)dt 



t 

The variance is calculated as follows: 

t+5T t+5T 

2 {{w-{wm 1 f r {limt")) 

^W- -77^ / dt I dt 

t+5T t+ST 



(VF)2 {5T)^ J J {ly 



^ J J dt" 



(I E{t') |2| l2\ 



J J {\Eit')\^){\Eit")\^) 

ST 5T 

du / dvg2{u -v)-l . (42) 







Substituting (|0|), and (|3|) into (||) we obtain the following expression for 
the variance of the energy distribution detected during time interval 6T for the 
Gaussian monochromator line profile: 



ST ST 



CTw = j^jy I du I dv\ gi{u-v) 





ST 



^ I (5T-T)Ui(r)|2dr 



-erf(5r(Tm) TFT^ • (43) 



5Tara {STan 



It can be shown that the integrated intensity distribution is described rather well 
by the gamma probability density function ([2l| ) with parameter M equal to cr^. 
When 6T is less than the coherence time Tc (|3^), the parameter M tends to unity 
and the gamma distribution tends to the negative exponential distribution. In the 
opposite case, when 6T ^ Tq, we can write: 



STam ST ' 

and the gamma distribution tends to a Gaussian distribution. 



In some experiments, such as one in which a long probe pulse is being used to 
measure lifetime of internal conversion of electronic excitation, the photocount time 
dependence is exponential dependence. In this case the detector gating function is 
an exponential function and the variance is calculated as follows: 
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2,0 



Fig. 9. Plot of M = a^J' versus (t^t^, exact solution for Gaussian spectral profile 



and exponential dependence of gating function 

oo oo 



/dt / dt'exp(-t/ro)exp(-t7ro)(| E{t) ^ E{t 



l\ |2\ 



2 







/ dtexp(-t/ro)(| m |2) / dt'exp(-tVro)(| E{t') |2) 



1 . 



(44) 



Let the monochromator line profile be Gaussian. Substituting (pO|), (^Tj) and ( |35| ) 
into (44) we obtain an expression for the variance of the energy distribution: 



oo oo 

= / d»e,p(-„)/d.exp(-„)e.p[-,^J(„- „)= 



(45) 



There is no simple analytical expression for the variance in this case. In Fig. ^ we 
present the results of numerical calculations. 

Our discussion of integrated intensity fluctuations would be incomplete if we did not 
refer the reader to another method of calculation. The sample dynamics is in the 
femtosecond range and (/) = const, is a good assumption. Clearly, we can assume 
that the synchrotron light in question is adequately modelled as an ergodic and 
hence stationary random process. Surely the preferred way to solve the problem of 
integrated intensity fluctuations must be stationary random process way. Although 
the origins of the above viewpoint are quite clear and understandable, the conclu- 
sions reached regarding the relative merits of stationary and nonstationary analysis 
are very greatly in error, for several important reasons. 

The second method of calculation has the advantage of conceptual simplicity but 
the disadvantage of being not quite as general as the first method. In defence of 
our nonstationary analysis we must say that models constructed for nonstationary 
analysis are inherently more general and flexible. Indeed, they invariably contain 
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the model of a stationary random process as a special case. In the normal course of 
events, a physicist first encounters statistical optics in an entirely stationary random 
process framework. In this case correlation functions in the frequency domain are 
represented by delta functions. The physicists emerging from such an introductory 
course may feel confident that they have grasped the basic statistical optics concepts 
and are ready to find a precise answer to almost any problem that comes their way. 
Nevertheless, there is no difficulty in compiling a long list of examples for which an 
analysis of nonstationary random processes is required. 

The nonstationary random process approach is needed somewhat more complex 
than the stationary random process approach, for it requires knowledge of the sta- 
tistical properties of synchrotron radiation not only in the time domain but also in 
the frequency domain. In the long run, however, the nonstationary random process 
model is far more powerful and useful than the stationary random process model in 
solving physical problems of practical interest. 



4 Photocount fluctuations in the detector of synchrotron radiation 

Almost all devices that measure the intensity of a light beam depend for their op- 
eration on the absorption of a portion of the beam, whose energy is converted to 
a detectable form. The photographic plate, the photomultiplier tube, and ioniza- 
tion and detection of photoelectrons in pump-probe experiments all fall into this 
category. The present section therefore considers the absorption of light in more 
detail in order to understand the nature of the intensity measurement process. It 
is convenient to base the discussion on a practical intensity-measuring device, and 
we select the photodetector. The photodetector depends for its operation on the 
photoelectric effect. Since the vast majority of pump-probe detection problems do 
indeed rest on the photoelectric effect, there is relatively little loss of generality by 
making this assumption. 

Consider a storage ring synchrotron light source. The light from this source falls on 
the photosurface, and we wish to determine the statistical distribution of the number 
of photoevents observed in any radiation pulse. Any measurement of light will be 
accompanied by certain unavoidable fluctuations. We attribute these fluctuations to 
quantum effects; that is, light can be absorbed only in quanta. We shall deal with the 
so-called semiclassical theory of photodetection. Such an approach has the benefit 
of being comparatively simple in terms of the mathematical background required, as 
well as allowing a greater tisc of physical intuition. The distinguishing characteristic 
of this formalism is the fact that electromagnetic fields are treated in a completely 
classical manner until they interact with the atoms of the photosensitive material 
on which they are incident. Thus there is no necessity to deal with quantization 
of the electromagnetic field; only the interaction of the classical field and matter is 
quantized. Fortunately it has been shown that the predictions of the semiclassical 
theory are in complete agreement wit the predictions of the more rigorous quantum 
mechanical approach for all detection problems involving the photoelectric effect 
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We assume that the radiation reaching the detector has full transverse coherence 
and that the statistical properties of the radiation follow the laws described in the 
previous section. It has been shown above that the energy, W, in the radiation pulse 
is unpredictable. Thus, we can predict the probability density p{W) only. In this 
case the probability of detection of K photons is given by Mandel's semiclassical 
formula 

P{K) = J eM-aW)P{W) d W , (46) 



where a = ij/hujQ and r/ is the quantum efficiency of the photodetector. It should 
be noted that, when conditioned by knowledge of energy W, the number of counts 
K is a Poisson variable with mean aW. Using formula ( ^6| ) we get the expression 
for the mean and for the variance of the value of K: 

{K)=a{W), 4 = ^^^^!^^ = ^ + ^^, (47) 

where a'^/ = 1/M is given by the expression (|2^). The expression for photocount 
fluctuations contains two terms. The first term corresponds to the photon "shot 
noise" and its origin is in the Poisson distribution. The second term corresponds to 
the classical fluctuations of the energy in the radiation bunch and takes its origin 
from the shot noise in the electron bunch. The ratio of the classical variance to the 
"photon shot noise" variance is named as a photocount degeneracy parameter 6c- 

The probability density of the integrated synchrotron radiation intensity, p(W), is 
the gamma distribution. Substituting (|26| ) into Mandel's formula ( ^6t) and perform- 
ing integration we come to the negative binomial distribution 

^ ^ r{K + i)r(M) \ {K)J \ M J ^ ^ 

When the monochromator has a narrow linewidth, AtOmTQ <C 1, parameter M tends 
to the unity, and the negative binomial distribution transforms to the Bose distri- 
bution: 

lim P{K) - ^ ' 



M^l ^ ' (l + (i^))A'+l • 

The negative binomial distribution tends to the gamma distribution at large values 
of the count degeneracy parameter 5c- In particular, the Bose distribution tends to 
the negative exponential distribution (5c — (K) in this case): 

(K)^ I f K \ 
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In the opposite case, at 5c — > 0, the negative binomial distribution ( ^9[ ) transforms 
to the Poisson distribution: 



K 



hm 

5c ^0 



Kl 



exp(-W) . 



5 Characteristics of undulator radiation 

For the proposed concept of time-resolved experiments based on the statistical prop- 
erties, the critical parameter determining the performance is the count degeneracy 
parameter, or the average number of photoevents produced in a single coherence in- 
terval of the incident radiation. The important role of this parameter is emphasized 
by giving it a name of its own. Because the count degeneracy parameter is propor- 
tional to {K) , it is also proportional to the quantum efficiency of the photosurface. 
Sometimes it is useful to remove this dependence on the characteristics of the par- 
ticular detector that may be present, and to deal with a degeneracy parameter that 
is a property of the synchrotron light itself. We thus define the wave degeneracy 
parameter as (5w = ^c/^- This new degeneracy parameter may be considered to 
be the count degeneracy parameter that would be obtained with an ideal detector 
having a quantum efficiency of unity. The wave degeneracy parameter is simply the 
average number of photons per mode. 

In this section we are particularly concerned about the wave degeneracy parameter. 
From this section we will see that the degeneracy parameter is (5w — A^-Bpeak/(4c), 
where 5peak is the peak spectral brightness of the radiation source. Third-generation 
synchrotron light sources achieve high spectral brightness through very low-emittance 
electron beams and by inserting into the electron beam a device known as an undu- 
lator. The undulator causes the electron beam to undulate or bend back and forth, 
and this produces very bright and coherent light. 



5.1 Number of photons radiated in the central radiation cone 

Let us introduce the basic features of undulator radiation. The undulator equation 



tells us the frequency of radiation as a function of undulator period Aw = 2Tik^, 
undulator parameter i^w? electron energy moc^'j, and polar angle of observation 6. 
Note that for radiation within the cone of half angle 



1 -1 
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the relative spectral FWHM bandwidth is ^ujjuj = O.SS/A'w, where is the num- 
ber of undulator periods. The spectral and angular density of the radiation energy 
emitted by a single electron during the undulator pass is given by the expression 
(at zero angle): 

^£ _ e^N^^'^Al^Kl sin2[7riVw(cj - a;o)M] 



du;dJ^ 2c(l + i^2/2)2 [vriV„(u; - wo)M]2 

Here = 2'^'^k^l{\ + is the resonance frequency, ^jj = [Jo(Q) — JiiQ)]-, 

where J„ is the Bessel function of nth order, and Q = K^/(4 + 2K'^). Now we 
would like to calculate the energy radiated into the central cone. In the small-angle 
approximation the solid angle is equal tod^^ = ^d^d9'- Integration of spectral and 
angular density over oj and 99 gives us factors ujq/N^ and 27r, respectively. We also 
have to integrate over 6 from to ^cen- Thus, the energy radiated into the central 
cone by a single electron is given by 

- c(l + Kl/2) ■ 

Individual positions of electrons in the bunch are random, thus the radiated fields 
due to different electrons are uncorrelated, and the average energy radiated by the 
bunch is a simple sum of the radiated energy from the individual electrons. 

Beyond the natural broadening, due to the finite number of oscillations, further 
spectral broadening can be incurred with the passage of many electrons through the 
undulator in a bunch of finite size, divergence, and energy spread Q. If there is an 
electron energy (rms) spread within the bunch, A7/7, there will be a corresponding 
photon energy spread given by AE'ph/E'ph = 2A7/7, where the factor of two is 
due to the square relationship between photon energy and electron energy. A more 
significant effect is that due to angular distribution within the bunch. As a result, 
some electrons traverse the undulator not along or parallel to the z-axis, but at 
a small angle Ox,y These electrons undergo the same number A^^ of oscillations, 
but experience a somewhat longer period A^/ co'8{6x,y)- If there is an rms angular 
divergence a' within the bunch, there will be a corresponding photon energy spread 
given by AE^ph/^Eph = 7^(o"')^/(l + In what follows we use the following 

assumption: 

A7 1 , „9 1 + Kll2 , , 

When these conditions are satisfied, the energy spread and angular divergence cause 
a spectral broadening less than and the central cone will be rather well defined 

in terms of both its angular definition and spectrum. 

To complete the calculation of the undulator characteristics an expression is needed 
for the photon flux within the central cone. We divide the radiated energy by the 
energy per photon hujo, and obtain 

dt I + Kl/2 ' 
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where a ~ 1/137 is fine-structure constant, N is the number of electrons in a bunch, 
/ is the bunch repetition rate. 



5.2 Spectral brightness of undulator radiation 



The quality of the radiation source is described usuahy by the spectral brightness 
defined as the density of photons in the six-dimensional phase space volume : 

where ax^y and {cTx^y)' are the (horizontal and vertical) photon beam Gaussian 
{l/^/e) radius and Gaussian half angle respectively. In the case of finite particle 
beam emittance the photon brightness is reduced. The amount of reduction, how- 
ever, depends on the matching to the photon beam. The particle beam parameters 
are 

i^e)x,y ~ {^e)x,yl3x,y , {'^e)x,y ~ {(^e)x,y/ Px,y ; 

where Px,y is the betatron function in the undulator, €x,y is the horizontal and 
vertical electron beam emittance, respectively. To modify the expression for spectral 
brightness for the case of undulator radiation, we can use the previously calculated 
photon flux in the central radiation cone, d-ZVcen/ dt, which was defined as having 
a relative spectral bandwidth (BW) of Auj/uq = 1/N„ (under condition (^)) and 
radiation cone of half angle ^cen- Taking into account the contribution of diffraction 
effects, we may wright the following expressions for the size and angular divergence 
of the photon beam: 

i^tot)l,y = {(^e)l,y + , {(^tot)l,y = {(^c)l,y + (^d)^ ) 

where = ^J XL^ / {%tt'^) is the diffraction limited radiation beam size, a'^ = 
\/ A/ (2Lw) is the diffraction limited radiation beam divergence, and is the un- 
dulator length 0. It is frequently used in the synchrotron radiation community to 
express the spectral brightness in terms of relative spectral bandwidth of 10~^. The 
expression for B in terms of photon flux within the central cone is Q : 

(dA^ccn/dt)(iVw/1000) photons/s 

ti — 



4vr2 (cTtot ) X (c^tot ) y (cr^ot ) x {(^[ot ) y mm2mrad^ (0.1 %B W) 



5. 3 Spatial and spectral filtering of undulator radiation 



Our previous discussion assumed that the synchrotron light striking the photosur- 
face was completely coherent in a spatial sense. In such a case the number of degrees 
of freedom M is determined strictly by temporal effects. When the wave is not spa- 
tially coherent, its spatial structure can affect the number of degrees of freedom; at 
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Fig. 10. Undulator radiation with a pinhole spatial filter 

any given time, different parts of the photosensitive surface may experience different 
levels of incident intensity. One of the way to overcome this problem is based on 
the idea to use a spatial filter. This possibility to create spatially coherent radia- 
tion is important for many experiments specifically for our scheme for time-resolved 
experiments and we will discuss in more detail the conditions for the synchrotron 
light source to emit such radiation. 



As an example of pinhole spatial filtering, Fig. 10 illustrates how the technique 
is used to obtain spatially coherent radiation from a periodic undulator The 
limiting condition of spatially coherent radiation is a space-angle product d ■ 6 = 
X/{2tt), where d is a Gaussian l/y^ diameter and 6 is the Gaussian half angle |§. 
In general the phase space volume of photons in the central radiation cone is larger 
than the limiting condition required for spatial coherence. Thus, for experiments 
that require spatial coherence, a pinhole and angular acceptance aperture are to be 
introduced. This pinhole spatial filter is used to narrow, or filter, the phase space 
of transmitted radiation. Filtering to d-9 = A/(27r) requires the use of both a small 
pinhole (d), and some limitation on 9, such that the product is equal to A/(27r). 

The radiation is transversely coherent when 

{cre)x,y(.0-'^)xy < A/(47r) , /3 ~ L„ . 

Under these conditions, one has (cr^.y)^ ^ ('^d)^ and {cr'x^y)^ ^ ('''d)^ ^'^'^ 
spectral-angular dependence of the radiation emitted by an electron beam can be 
approximated by the spectral-angular dependence of the radiation emitted by a 
single electron. This limit corresponds to the maximum spectral brightness 

^ ^ 4(d^en/ d i)(iV„/1000) photons/s 



A2 mm2mrad2(0.1%BW) 
When the electron beam emittance is large, 
> A/(47r) , 
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Fig. 11. Young's interferometer, for discussion of spatial coherence 

the radiation is partiahy coherent. One can show that the influence of the electron 
beam divergence on the properties of the undulator radiation can be neglected when 
the beta-function is large enough. In the region of parameters when 

a downstream diaphragm of aperture d = 2a ^ is used for selection of the transversely 
coherent fraction of undulator radiation. Finally, in this limit, the flux of transversely 
coherent photons into the bandwidth Aw/w = 0.1% can be estimated simply as 

diVph 
dt ^ 4 ■ 

These remarks may be useful in estimating the conditions under which a spatially 
coherent plane wavefront may be obtained from a noncoherent source. 



5.^ Discussion of spatial coherence 



The spatial coherence in a light beam generally has to do with the coherence between 
two points in the field illuminated by light source. The meaning of spatial coherence 
can best be understood with the help of Young's two-pinhole experiment (see Fig. 



11). In its elementary sense, the degree of coherence between the two points simply 
describes the contrast of the interference fringes that are obtained when the two 
points are taken as secondary sources. Let a source S illuminate the two pinholes 
Si and S2, as shown in Fig. |ll[ The source is perfectly noncoherent. That is to say, 
no interference fringes can be obtained by placing two pinholes in the plane of the 
source. It was shown, however, that if the two slits are placed far enough away from 
the noncoherent source, interference fringes of good contrast can be obtained (see 



Fig. 12). It is sometimes said that the spatial coherence in light beams increases 
with distance "by mere propagation". It would be nice to find an explanation which 
is elementary in the sense that we can see what is happening physically. 

Suppose that a quasimonochromatic wave is incident on an aperture in an opaque 
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Fig. 12. Interference-fringe amplitude, in Young's interferometer, for three values of 
noncoherent source-pinhole screen distance 
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Fig. 13. Geometry for propagation of spatial coherence 



screen, as illustrated in Fig. 13. In general, this wave may be partially coherent. 
The detailed structure of an optical wave undergoes changes at the wave prop- 
agates through space. In a similar fashion, the detailed structure of the spatial 
coherence undergoes changes, and in this sense the transverse coherence function is 
said to propagate. Knowing the spatial coherence on the aperture, we wish to find 
the spatial coherence on the observing screen at distance z beyond the aperture. 
Synchrotron radiation is a stochastic object and for any synchrotron light beam 
there exist some characteristic linear dimension. Arc, which determines the scale of 
spatially random fluctuations. Fig. |l^ illustrates the type of spiky speckle pattern 
on an aperture in an opaque screen. When Arc ^ the radiation beyond the aper- 



ture is partially coherent. This case is shown in Fig. 13. Here Arc may be estimated 
as the typical linear dimension of speckles. 



First we wish to calculate the (instantaneous) intensity distribution observed across 
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Fig. 14. Reciprocal width relations of Fourier transform pairs 



a parallel plane at distance z beyond the aperture. The observed intensity distri- 
bution can be found from a two-dimensional Fourier transform of the field. The 
radiation field across the aperture may be presented as a superposition of plane 
waves, all with the same wavenumber k = ujo/c. The value of k_i_/k gives the sine 
of the angle between the z axis and the direction of propagation of the plane wave. 
In the paraxial approximation k±/k = sinO ~ 0. If the radiation beyond the aper- 
ture is partially coherent, a spiky angular spectrum is expected. The nature of the 
spikes in the angular spectrum is easily described in Fourier-transform notations. 
We can expect that the typical width of the angular spectrum envelope should be 
of the order of (/cArc)~^. Also an angular spectrum of the source having transverse 
size d should contain spikes with a typical width of about (kd)^^, a consequence 



of the reciprocal width relations of Fourier transform pairs (see Fig. 14). It is the 
source linear dimension d that determines the coherent area of the observed wave 
z/{kd), but in addition the coherence linear dimension Arc of the source influences 
the distribution of average intensity over the observing screen with typical width 
z/{kArc)- Thus, if the screen is placed far enough away from the noncoherent source, 
z ^ dAvc/X, a coherence area of large linear dimension can be obtained. 



5. 5 Grating monochromator considered as a spatial filter 



In order to use the radiation from the source it first has to be filtered and monochro- 
matized. Finally, to be complete, we should remark that in the experimental ar- 
rangement shown in Fig. ^ there is still one other effect which gives spatial filtering. 
A monochromator is just a spatial filter working in dispersive direction only. We 
would like to show next how it is that a monochromator can act like a spatial filter. 



Figure 15 shows the grating geometry under consideration. The finite size of the 
entrance and exit slits in the dispersive direction sets a limit to the achievable res- 
olution. The slit-width-limited resolution can be obtained directly from the grating 
equation. In particular, the entrance slit-width-limited resolution is given by Q: 

AA = zdcos{a)/{nm) , 
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Fig. 15. niustration of grating monochromator. The finite size of the entrance sHt 
in the dispersive direction sets a hmit to the achievable resolution 



where a is the angle of incidence, A'g is the number of lines on the grating, D is 
the width of grating, n = Ng/D is the line density, m is the order of the diffraction 
pattern, d is the width of entrance slit, z is the entrance slit-grating distance. We 
see that the resolution is equal to 

A zmNgX 
AX Ddcos{a) 

The coefficient {mNg)~^ is just the minimal uncertainty in the wavelength that 
can be measured with the given grating. The coefficient of proportionality between 
A/AA and mNg is zX/[Ddcos{a)]. What does all this mean? This is an example in 
which transverse coherence plays an important role. The ability to spectrally filter 
radiation by a monochromator requires well-defined phase and amplitude variations 
of the fields across the grating. So in order that we shall have a sharp line in 
our spectrum corresponding to a definite wavelength with an uncertainty given by 
{mNg)~^, we have to have a wave transverse coherence area of at least the size 
Z)cos(a). If the transverse coherence area is too small we are not using the entire 
grating. Suppose we have an entrance slit with aperture d and send a beam of 
radiation at the grating with an aperture -Dcos(a). The radiation reaching the 
grating is transversely coherent when the following condition for the space-angle 
product is fulfilled: 

d-e = Ddcos{a)/{2z) ~ X/{2tt) . 

We expect that only for pure transverse coherence, the grating will work right and 
the resolution A/AA and mNg will be nearly equal. Importantly, in our scheme that 
incorporates an external pinhole spatial filter (see Fig. |^, the monochromator does 
not require an entrance slit. 
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5. 6 Wave degeneracy parameter of undulator radiation 



Physically the wave degeneracy parameter describes the average number of photons 
which can interfere, or, according to quantum theory, the number of photons in one 
quantum state (one "mode"). We shall deal with two different approaches to the 
calculation of the wave degeneracy parameter. First we consider what can reasonably 
called "fluctuation experiments in the frequency domain" (see subsection 3.6). For 
the frequency domain we use the notion of the wave degeneracy parameter 5w 
which is equal to the average number of transversely coherent photons radiated 
by the electron bunch inside the spectral interval of coherence Awc- Taking into 
account that the "emittance" of the diffraction limited photon beam is equal to 

mm.{ex,y) = min(27rcr^,j^c7^ j^) = , 

we calculate the average number of spatially coherent photons radiated within one 
pulse into the spectral interval of AA/A: 

A^AA 



A AA f 

(AiVph)coh = XX i ' 



where B{t) is the instantaneous spectral brightness of the synchrotron radiation. 
The value of AA/A is connected with the interval of the spectral coherence by the 
relation: AA/A = Awc/wo- Using (p^), (p^) and Parseval's theorem we obtain: 

oo oo oo 

^ = - / I 5i(Aa;) |2 d(Aa;) = - / | F(Au;) ^ d(Au;) = - / F\t)dt. 

— oo — oo — oo 

Using the fact that the instantaneous spectral brightness B{t) and the instantaneous 
value of the beam current {I{t)) are related by B{t) oc {I{t)), we can equivalently 
write B(t) = const, x F{t). It is convenient to refer to the peak value of spectral 
brightness 

max(i?) = const, x max(F) = -Bpeak • 

By convention, we represent the time integral of the instantaneous spectral bright- 
ness in the form 

oo oo 

/B(*)<it = const.x/F(t)dt = ^. 

— OO — oo 

Note that the dependence of the factor 

oo 
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on the exact shape of the bunch is rather weak. The results are A = 1 for the 
rectangular pulse-shape and A = l/\/2 for the Gaussian pulse-shape. Thus ^ ~ 1 
is a reasonable approximation. Finally, the degeneracy parameter can be estimated 
simply as: 

<5w ^ . (52) 

What about the other kinds of fluctuation experiments, for example, integrated 
intensity fluctuations? For the time domain we use the notion of the wave degeneracy 
parameter 5w which is equal to the average number of transversely coherent photons 
radiated by the electron bunch during the coherence time Tc- 

(AiVph)coh = (''"c-Bpcak) • 

Finally, the product TcAX/X is given by the expression 
TcAA/A = A'\/c . 

The dependence of the factor A' on the exact shape of the spectral distribution 
is rather weak. Let us consider the specific case of a Gaussian line profile. The 
explicit expression for the coherence time is = -v/vr/<^m- The FWHM bandwidth 
would be AA/A = 2.36((Tni/'^) giving, for A' = 0.7. Thus A' ~ 1 is a reasonable 
approximation and the wave degeneracy parameter is given by the expression 



(52). The frequency domain calculations and time domain calculations give similar 



results. 

Let us present a specific numerical example for the case of a third-generation syn- 
chrotron light source. The peak spectral brightness at a wavelength of 10 nm is 
equal to 



mm2mrad2(0.1%BW) 



Substituting this number into (^2]) we obtain that the wave degeneracy parameter is 
about (5w — 10^. Hence in the soft X-ray region of the spectrum we expect classically 
induced fluctuations of photocounts to have a far stronger effect than pure quantum 
shot noise fluctuations. On the other hand, in the hard X-ray region of the spectrum 
(A ~ 0.1 nm), a source brightness in excess of 

5peak - 1024 Photons/s 

P mm2mrad2(0.1%BW) 

is required to produce a wave degeneracy parameter greater than unity. Since third- 
generation sources have a peak brightness of only this value, we can conclude that in 
this region of the spectrum, the vast majority of third-generation synchrotron light 
sources encountered produce radiation with wave degeneracy parameters compara- 
ble with unity, and hence noise produced by the quantized nature of the radiation 
is comparable with the noise produced by classical fluctuations of the intensity. 
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A few additional comments are needed in closing this section. We have concen- 
trated our attention on the wave degeneracy parameter, which is a property of the 
synchrotron radiation falling on the photodetector. That photodetector invariably 
has a quantum efficiency that is less than unity. In addition, the efficiency of the 
monochromator is less than unity, too. Hence the count degeneracy parameter will 
be even smaller than the wave degeneracy parameter. When the integration time is 
large, tq 3> Tc, the parameter is given by: 

(5c = r?-Rm<5w ■ 

where is the monochromator throughput. In the opposite case, at tq <C Tc, the 
parameter 6c is about 

6c ^ r]Rm{^LOmTo)6w ■ 

In addition, it is possible that the sample or collecting optics may intercept only a 
fraction of one spatial mode from the source. In such a case the count degeneracy 
parameter may again be smaller than the wave degeneracy parameter, as a result 
of the incomplete capture of a spatial mode. Although the minimum value of M 
is unity, the reduction of the energy striking the detector surface must be taken 
into account. In such a case, the count degeneracy parameter must be reduced from 
the normal value by a factor of the ratio of the effective measurement area to the 
coherence area of the incident light. 



6 The signal-to-noise ratio associated with the output of the mea- 
surement system assumed for the pump-probe experiments 



The principle of operation of the proposed pump-probe scheme is based essentially 
on the classical intensity fluctuations of synchrotron radiation. It is shown above 
that stochastic fluctuations of the classical intensity can influence the statistical 
properties of the photoevents that are observed. Note in particular that the variance 
(T^ of K consists of two distinct terms, each of which has a physical interpretation: 

^K = ^^y^ = ^ + ^w (53) 

The first term is simply the variance of the counts that would be observed if the 
classical intensity were constant and the photocounts were purely Poisson. We re- 
fer to this contribution to the count fluctuations as "quantum noise". The second 
term, a^^, is clearly zero if there are no fluctuations of the classical intensity. For 
instance, in the case of laser light, this component would be identically zero, and 
the count variance would be simply that arising from Poisson-distributed counts. 
When synchrotron light is incident on the sample surface, the classical fluctuations 
are nonzero, and the variance of the photocounts is larger than that expected for 
a Poisson distribution by an amount that is proportional to the variance of the in- 
tegrated intensity. This extra component of variance of the counts is often referred 
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to as "excess noise," meaning that it is above and beyond that expected for pure 
Poisson fluctuations |Q. 

Our initial goal was to find expressions for the classical variance cr^ of the integrated 
intensity. Also of major interest is the signal-to-noise ratio, associated with the 
variance, which provides us with an indication of the magnitude of the fluctuations 
of classical variance relative to the expectation value. The general approach to 
calculating the output signal-to-noise ratio of the proposed device will be as follows. 
Following subtraction of the {K) associated with the quantum noise (see (|53|)), the 
signal-to-noise ratio takes the formll 



where cr^ is calculated using (^), N^hot is the number of independent measurements 
averaged in the accumulator (total number of shots). The only requirement for the 
accuracy of this procedure is that the count fluctuations be uncorrelated from shot 
to shot, a property that does hold in our case. 



A fully general study of the signal-to-noise ratio would be nontrivial. The difficulty 
arises in simultaneously including the effects on that noise of both the classically 
induced fluctuations and Poisson fluctuations of the counts. The full analysis, in- 
cluding both of these effects, is a difficult analytical problem. The calculation can be 
performed without great difficulty in two limiting cases, namely, the cases of degen- 
eracy parameter very large and very small compared with 1. When (5c ^ 1, there 
are many photoevents present in each coherence interval of the wave. The result 
is a "bunching" of the photoevents by the classical intensity fluctuations, and an 
increase of the variance of the counts to the point where the classically induced fluc- 
tuations are far stronger than the quantum noise variations. When the degeneracy 
parameter is much larger than 1, the expression (|5^ ) can be replaced by 

S VN^.{(K - (K)r-) 

N - (A-))«) - UK - {K)W 

Our consideration is restricted to the two simplest cases from the analytical point 
of view, namely, the cases of integration time tq very short and very long com- 
pared with the coherence time. For short integration time, the value of the incident 
intensity I[t) is approximately constant over the entire counting interval. As a con- 
sequence, the integrated intensity W is distributed in accordance with the negative 
exponential distribution (pO|). It follows that the (-ftT") = ((aVF)") can be expressed 



^ A complete analysis of the finite averaging should include the uncertainties asso- 
ciated with the estimates of all average quantities. For the purpose of simplicity we 
neglect the uncertainty of {K). An assumption that this quantity is known much 
more accurately than {{K — {K))'^) is justified in our case, actually several or many 
different coherence times must be explored, and {K) is thus observed over many 
more shots than is {(K — {K))'^) for any one coherence time. 
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in terms of the mean as 
= n! . 

Using the above expression, we find that the signal-to-noise ratio is given by the 
expression: 

S/N ~ ^A^shot/S for 5c > 1, > tq . (56) 

This calculation shows that the fluctuations of the {{K — {K))"^) is about -y/S ~ 3 
times larger than the fluctuations of the mean value {K) . 

At the opposite extreme, with an integration time much longer than the coherence 
time, the integrated intensity, W, is distributed in accordance with the Gaussian 
distribution. It follows that 

S/N ~ for 5c » 1, tq > . (57) 

We know that for the case when 5c <^ 1 the fluctuations of photocounts are strongly 
dominated by pure quantum noise. Thus the probability of obtaining K counts is, 
to a good approximation, given by the Poisson distribution. We cannot neglect the 
classically induced fluctuations of the counts when we calculate the signal compo- 
nent of the output, but we can neglect them when we calculate the noise, simply 
because their contribution to the noise is so small. For the final approximation we 
note that if the integration time tq is very short compared with the coherence time, 
the expression for the photocounts degeneracy parameter is equal to the average 
number of detected photoelectrons per pulse 5c — {K). Hence, when 5c ^ 1 and 
To <C Tc the following approximation is valid : {K) ^ 1. As a reminder, recall that 
the factorial moments of the Poisson distribution are given by 

{K{K -1)---{K -n + l)) = {K)"" . 

It follows that 

{{K - {K)f) - {{K - {K)f)^ ~ {K) for {K) « 1 . 

Substituting this result into ( [5^ ) we find 

S/N ~ 53/2 ViV^ for <5c « 1 , TO « Tc . 

Note that this expression is independent of the integration time tq. Therefore, the 
signal-to- noise ratio is not improved by using a longer " characteristic time" of sam- 
ple dynamics. 



7 Pump-probe technique based on a correlation principle 

A more complex instrument can be used to characterize the shape of the signal 
function. The technique we are to describe here is to use two spatially separated 
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Fig. 16. Scheme for pump-probe experiments employing correlator of count fluctu- 
ations 




Fig. 17. Correlation apparatus for analysis of sample dynamics 
samples and record correlation of the count fluctuations for each value of the delay r 



between two pump laser pulses. The new configuration is sketched in Fig. 16. In this 
case the transversely coherent probe X-ray beam is splitting and incidents on two 
spatially separated samples. The detected photoevents are than correlated, and the 
sample dynamics is determined from that correlation. Each detector is connected 
with a separate electronic counter. The version of apparatus for measurement of 



average count product depicted in Fig. 17. During one radiation pulse, each of 



two detectors registers Ki and K2 number of photocounts, respectively. After each 
shot, an electronic scheme multiplies these numbers and passes this product to an 
averaging accumulator, where it is added to the previously stored sum of count 
products. Finally, the total sum is divided by the number of shots. This result, 
averaged (shot to shot) count product {K1K2), contains information about the 
shape of the signal function. Another advantage of the correlator measurement is 
the possibility to remove the monochromator between spatial filter and sample and 
thus to increase the count degeneracy parameter. 

Now our goal is to find the relationship between the expected value of the average 
count product and the sample dynamics. In addition, we wish to find the variance 
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of this quantity, so that the signal-to-noise ratio associated with the measurement 
can be determined and compared with the similar quantities found in the previous 
section. 



7. 1 The expected value of the count-fluctuation product and its relationship 
to the sample dynamics 



By the "count fluctuations" we mean explicitly the difference between the actual 
numbers of counts obtained in one shot at detector 1 and 2 and the expected values 
of these two number counts. Thus 



AKi = Ki-{Ki) , AK2 = K2- {K2) 



The averaging accumulator at the output of the system depicted in Fig. 17 in effect 
produces an estimate of the expected value of the product of the two count fluctu- 
ations. Thus we are interested in statistical properties of the quantity (Ai^'iA^C2)) 
and particular its mean and variance. 

Let us calculate the expected average value of {K1K2): 

00 00 

{K1K2) = E E KiK2P{Ki,K2) , 

/Ci=0 ^2=0 

where P{Ki, K2) is the joint probability distribution of Ki and K2. It follows from 
the basic properties of the conditional probabilities that 

00 00 

P{Ki,K2) = JdWi JdW2P{Ki,K2 I Wi,W2)PiWi,W2) , 


where P{Wi, W2) represents the joint probability distribution of integrated intensi- 
ties Wi and W2- The values Wi^2 can be interpreted as the portion of X-ray pulse 
energy converted to a detectable form - photoelectrons for example (see section 2.7). 
Since Ki and K2 are independent when conditioned by the integrated intensities 
Wi and respectively, we can write Q: 

P{Ki,K2 I Wi,W2) = P{Ki I Wi)P{K2 I W2) , 

Using Mandel's formula (p^ ) we can write: 



P{K I W) = ^— ^exp(-aW) . 
K'. 

Using these facts, the following expression for the average of the count product can 
now be written 



{K^K2) = V V KrK2 / / ^ J{ e^-aWi Y J\ expC-aT^a) 



Xi=0i^2=0 



Ki\ ^' K2I 
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xP{Wi,W2)dWidW2 



(58) 



At this point we interchange the orders of summation and integration in (|58| ) . Using 
the relationship 

^ ' exp(-aVF) = aW , 

K=0 

we can express the average of the count product in terms of the average of the 
classical integrated intensities at the two detectors Q, 

{K1K2) = a^{WiW2) . 

This result shows that much can be learnt about the classically induced fluctuations 
by studing the photocount correlation. So, the next problem is the description of 
the correlation of the classical integrated intensities 



Wi,2= j h,2{t)h,2{t)dt 



where f{t) is the sample gating function (signal function). Suppose that the pump 
pulse interacts with the first and second sample at t = and t = — r, respectively. 
In this case the gating functions are: fi{t) = f{t), f2{t) = fit + ''")• When the 
average intensities incident on the two samples are constant and equal, 

{h{t)) = {hit)) = (/)= const. , 

the expression for the correlation of count fluctuations reduces to the useful form! 



{AK1AK2) {AW1AW2) 



{Ki){K2) {Wi){W2) 



1 



'ch 



00 00 



J dt J dt'fit)fit' + T) 



-00 —00 
00 00 



mm) 



^ [ dt [ dt'fit)fit' + T)\g,it-t') |2 

^ch 



(59) 



-00 —00 



where r^h denotes: 

00 

Tch= J fit)dt. 
— 00 

^ This case is particularly interesting, because sample dynamics is in the femtosec- 
ond range, and (/) = const, is a good assumption. 
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signal function autocorrelation function 

versus time versus delay 

f(t) A^(T) 




t T 



Fig. 18. Some typical signal versus time and corresponding autocorrelation Ac{t) 

Then we remember that the behavior of | gi{t — t') ^ for Tch ^ Tc approaches the 
behavior of the delta function (here Tch is the "characteristic time" of the sample 
dynamics). As a result we obtain 



{AK1AK2) 



{Ki){K2) 
One can also write (IBOD as 



00 

^ / /(t)/(t + T)dt for Tch>rc. (60) 

— CO 



00 



{AKIAK2) 



a 



I f{t)f{t + T)dt 

2 —CO 



7/^(t)dt ' 

—00 

where is the variance of the integrated intensity distribution ( for example see 
expression (^5|)). In the limit of a large photocount degeneracy parameter, 6c ^ 1, 
we have ~ cr^. 

The expression 
00 

Mr)= I f{t)f{t + T)dt 
—00 

is called an autocorrelation. One immediately recognizes the physical meaning of 
the autocorrelation functions. The Fourier transform of the autocorrelation is Ac{ijo)., 
related to the Fourier transform of the signal by: Ac{uj) =| f{oj) p. An autocorrela- 
tion is always a symmetric function. The Fourier transform of the autocorrelation 
is a real function, consistent with a symmetric function in the time domain. The 
question then naturally arises as to exactly what information about f{t) can be 
derived from the measurement of cross-correlation (see Fig. [l^ ). One can see that 
the proposed correlation technique provides the possibility to measure the modulus 
of the Fourier transform of the signal function, while information about its phase is 
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missing. To provide reconstruction of the signal function from the autocorrelation 
function a special iterative technique should be used. Clearly, without loss of gen- 
erality, we can assume that the signal function describing the sample dynamics is 
zero for t < and nonzero for t > 0. Additionally we know that at t = the signal 
function is maximum. These relationships help us with the task at hand, namely 
determining the phase of /(w) from the knowledge of its modulus in the majority of 
cases involving functions f{t) of considerable complexity. Nonetheless, it is possible 
to find cases in which ambiguities exist. 



7.2 The signal-to-noise ratio associated with the output of the correlation 
apparatus 



The signal-to-noise ratio for the proposed device (see Fig. [T^ ) can be written in the 
form: 



where N^hot is the number of independent measurements averaged in the accumu- 
lator. Let us consider the case when the degeneracy parameter is much larger than 
unity, 6c ^ 1- The fluctuations of the photocounts are defined mainly by the classical 



noise in this case and ( plD can be reduced to 



Remembering that the integrated intensity is distributed in accordance with the 
Gaussian distribution one can derive that 



2' 1 I ^2 



{AWiAW2) _ 2 , . ^ 

{W,){W2) - ^^^'^^^^ ' 



2 



(^W,)^W2)' = 1 + 2<T^ + ^'^wMr) + 2a^Al{T) + , (63) 



where A^{t) denote 

00 

/ f{t)f{t + T)dt 

Mr) = 

/ Pit)dt 
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Combining the above expressions, we find that the signal-to-noise ratio is given by 
the expression: 



S_ ^ VN~^tAniT) 



for 



> 1 



(64) 



Let us analyse the signal-to-noise ratio taking into account the quantum effects. 
When calculating the signal component (numerator of (|6l|)), we can take into ac- 
count classically induced fluctuations of the counts only. This can be done due to the 
fact that the quantum fluctuations of the counts at the outputs of two detectors are 
statistically independent. General calculations of the noise fluctuations associated 
with the output of the correlation apparatus (denominator of (|6l])) should include 
both the classical and the quantum effects. Taking into account these considerations 
we can derive the following relations: 



{aWi)^{aW2) 



{Ki)HK2)^ 
{KIK2) {{a^Wl + aWi)aW2 



(Kf^) {a^Wl2 + aWi,2) 



1,2 



(65) 



We illustrate the technique for calculations of the above expressions for the value 
of {KfKl): 



00 00 



{k!kI)= Y: E K!KlPiK„K2) 

Ki=0K2=0 



dWiJ dW2PiWi,W2) X 



Ki=0 
00 



00 

E K2 

K2=0 



2 iaW^2)^^ 



K2! 



exp(— q1^2) 



J d^^J dW2P{Wi, W2){a^Wl + aWi){a^Wi + aW2) 



(66) 



We will consider only the special case when the degeneracy parameter is much less 
than 1. Using the approximation (5c ^ 1 we obtain (5c) 1,2 — (-^^1,2)- Under such 
condition one would have, from (|65|), the variance 

{{AKiAK2f) - {AKiAK2f ~ {K1K2) , 
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and the signal-to-noise ratio would be 



|: = ^N^tScAnir) for ,5c « 1 . (67) 



The analysis of the obtained asymptotic expressions (^4|) and ( p7\ ) for the signal-to- 
noise ratio allows one to make the following conclusions: 

(i) The signal-to-noise ratio depends on the degeneracy parameter only when 5c is 
much less that 1. In this case the number of independent measurements for a given 
signal-to- noise ratio is proportional to (l/5c)^- 

(ii) The signal-to- noise ratio is proportional to A^-^^t). To hold the signal-to-noise 
ratio constant, the number of independent measurements must be proportional to 



7.3 Other autocorrelation techniques 



There are other techniques to provide the possibility to measure the integrated 
intensity autocorrelation. They differ from each other in the beam splitters, the 
delay line and the detection system. We will now describe two new approaches. An 
alternative technique to generate An{T) is to use one pump beam and two probe 
beams. In the setup shown in Fig. |l^ the pump pulse perturbs the sample at time 
t = 0. A long X-ray pulse passes a beam splitter, generating two pulses travelling 
different paths. After the two paths are recombined, the two pulses impinge on the 
sample at a time difference r which can be varied. The photocounts are recorded 
by a detector and the variance cj^ is studied as a function of r. 

Now we describe the basic principles of autocorrelation measurement with this ar- 
rangement. In the limit of a large photocount degeneracy parameter the expression 
for the variance of counts reduces to the form ((/i) = {I2) = const.): 



( {[hit) + hit + r)][h{t') + hit' + r)]) 



\ I dt f dt'fit)f{t')[- 



{{h+h)? 

It I I dt'fit)fit')[\giit-t')\^ + \giit-t' + T)\^] . (68) 



-00 —00 
00 00 



'CO —00 



Remembering that the behavior of | giit — t') p for Tch ^ Tc approaches the behavior 
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of the delta function we obtain El 
1 



^ch 

— oo 



[l + ^n(r)] for Tch>rc. (69) 



One can also write 



^wW = ^^w(0)[l+^n(T)] , 



where o"^(0) is the variance of the integrated intensity distribution for the single 
probe beam: 



oo 

— oo 



At T = 0, the peak value of the function (t'^{t) is ^^^^(O). For large delays compared 
to the "characteristic time" of the sample dynamics, the cross product term vanishes, 
leaving a background of cj|^(oo) = o"^(0) /2. The variance as a function of delay time 
has thus a peak to background of 2 to 1. The measurement leading to o^^ij) we can 
refer to as the integrated intensities autocorrelation with background, as opposed to 
the background free autocorrelation described in the previous subsection (see Fig. 



The most important and critical area of R&D for time-resolved experiments utilizing 
the autocorrelator will be the development of X-ray beam splitters and delay lines. 
Considerable work will need to be undertaken to make them routinely used optical 
elements. 

The technique to use two probe X-ray beams that are generated by the same electron 
bunch but from two insertion devices (again with one beam delayed with respect 
to the other) may be a more promising approach. The concept is very interesting 



in view of the critical synchronization of probe beams (see Fig. 2C). Using this 
technique precisely controlled temporal resolution down to several femtoseconds 
should be achievable. To delay the electron beam, a small magnetic chicane has to 
be applied. The trajectory of the electron beam in the chicane has the shape of an 
isosceles triangle with the base equal to L. The angle adjacent to the base, 0, is 
considered to be small. The optimal parameters of the chicane can be calculated in 
the following way. The electron beam delay line needed for time resolved experiments 
has to satisfy two requirements. First, the electron beam extra path length, L0^/2, 
must be long enough to reach a delay time of about 100 femtoseconds. A first 
question that arises is: why do we want to use 100 fs delay time? Here we direct 



^We can assume that r Tc, so that no interference is possible between waves 
of different probe beams. In this case the contribution to the light intensity is 
obtained by a suitable summation of the intensity contributions. Since we assume 
that Teh S> Tc, there is relatively little loss of generality by making this additional 
assumption. 
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delay 




Fig. 19. Scheme for autocorrelation pump-probe experiments employing two probe 
beams originating from the splitting of a single X-ray beam produced by an undu- 
lator (with one partial beam delayed relative to the other) 




magnetic chicane 

Fig. 20. Scheme for autocorrelation pump-probe experiments employing two probe 
beams that are generated by the same electron bunch but from two undulators 

the attention of the reader to the fact that we assume to use 5-10 fs visible pump 
pulses. On the other hand we assume that the pump pulse can be approximated 
by a delta function. One can never hope to achieve a visible laser pulse duration 
ten times shorter than the present. This limitation prevents the application of this 
technique to the timescale shorter than 100 fs. 

Second, the effect of suppressing the beam density modulation due to the presence 
of the energy spread in the electron beam should be avoided in order to preserve 
the correlation of the radiation field within the delay time. The energy spread sup- 
presses significantly the microbunching when ujqLO'^ A'y / {c^) is larger than or equal 
to unity. Parameters in our case are: A7/7 ~ 10~^, L6^ ~ 60/xm. The delay time 
and energy spread in the electron beam dictates the choice of synchrotron radiation 
wavelengths. These results demonstrate that this scheme is adequate for the VUV 
spectral range, but cannot be used to produce X-ray probe beams. The situation is 
quite different for the first scheme of autocorrelator (see Fig. ^). Since in this case 
the probe beams are produced by X-ray splitter and delay line, the energy spread 
limitation does not exist at all. 
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8 Concluding remcirks 

A general objective in the development of synchrotron radiation sources is to pro- 
duce radiation that is brighter than that from existing sources, or to produce ra- 
diation that comes in shorter pulses. Current machines produce bunch lengths in 
the 100-ps regime. Many experiments exploit this time dependence. Light-triggered 
time-resolved studies with 100 ps resolution are now frequently performed at syn- 
chrotron light sources. However, 100 picoseconds is longer than the time associ- 
ated with many interesting physical phenomena. The storage ring technology itself 
approaches its limits of performance with respect to minimal pulse duration. No 
practical means have been found to allow bunches of less than a ten picoseconds 
duration to be stored. A proposed in this paper schemes for fast time-resolved ex- 
periments could offer the possibility to overcome this limitation. It was shown that 
the femtosecond timescale associated with intensity fluctuations in synchrotron light 
pulses make them well suited for time-resolved studies. The new concept eliminates 
the need for ultrashort X-ray pulses. This development in synchrotron based pump- 
probe experiments allows us to investigate phenomena in the time range down to 
100 fs. It is the author's hope that the general principles of pump-probe techniques 
based on statistical properties of synchrotron light, presented in this paper, may help 
in stimulating further inquiry and result in success in many exciting applications, 
some of which can now only be barely envisioned. 
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